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VECTOR CALCULUS 


DIRECTIONAL DERIVATIVES AND GRADIENTS 


The partial derivatives 0z/dx and dz/dy are the rates of change of z = f(x, y) as the 
point (x, y) moves in the direction of the x-axis and the y-axis. We now consider the 
rate of change of z as the point (x, y) moves in other directions. 

Let P(a, b) be a point in the (x, y) plane and let 


U = cos ai + sin gj 


be a unit vector, « is the angle from the x-axis to U (see Figure 13.1.1). The line through 
P with direction vector U has the vector equation 


X=P+1U 
or in parametric form, 


() x=a+tcosa, y=b+tsina. 


Figure 13.1.1 The unit vector 


At t = 0 we have x =a and y = b. If we intersect the surface z = f(x, y) with the 
vertical plane through the line (Equation 1), we obtain the curve 


z=f(at+tcose,b+tsina) = Fit). 
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The slope dz/dt = F’(0) of this curve at t = 0 is called the slope or derivative of f in 
the U direction and is written f(a, b) (Figure 13.1.2). 


i directional derivative 


Figure 13.1.2 


Here is the precise definition. 


DEFINITION 


Given a function z = f(x,y) and a unit vector U = cosai + sinaj, the 
derivative of f in the U direction is the limit 
f(a + tcosa,b + tsina) — f(a, b) 


Sula, b) = lim P 
170 t 


Sula, b) is called a directional derivative of f at (a, b). 


The partial derivatives of f(x, y) are equal to the derivatives of f(x, y) in the 
i and j directions: 


Pad) tig LEO = KD) 


Ax>0 Ax 
ia J(a + tcos0,b + tsin0) — f(a, b) =a Fahy 
0 t 
Fabs itn /(a,b + Ay) — fla,b) 
Ay70 Ay 
flat tcos5.b + isin 5 — f(a, b) 
= lim _ = fila, 6). 


170 tf 


EXAMPLE 1_ Find the derivative of f(x, ¥) = xy + y? in the direction of the unit 
vector 
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sle+ Say + by) - T(x y) 


ful y) = lim 
170 t 


+ a [» + st)4 b + al — (xy + y’) 


= lim 
t-0 


Bra 


yt + e+ yt ae 


= fae 
1-0 


aii eg Broo 
> 


1 J3 
== — + 1]y. 
a+ [P+ i)y 


r 


There is an easier way to find the directional derivatives of f(x, y) using the 
partial derivatives. It is convenient to combine the partial derivatives into a vector 
called the gradient of /. 


DEFINITION 


The gradient of a function z = f(x, y), denoted by grad z or grad f, is defined by 


In functional notation, 


grad f = f(a, b)i + f,(a, b)j. 


Thus grad f is the vector valued function of two variables whose x and y 
components are the partial derivatives f, and f, (Figure 13.1.3). Sometimes the 
notation Vf or Vz is used for the gradient. 


Figure 13.1.3 


THEOREM 1 


Suppose z = f(x, y) is smooth at (a, b). Then for any unit vector U = cos ai + 
sin aj, the directional derivative fy(a, ne exists and 


z 
Sue, b) = U-grad f = 5 C08 + oy sin &. 
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PROOF Let U = cos ai + sin aj. Write x, y, and z as functions of f, 


x=a+tcosa, yeattsing, 
z=f(a+tcosa,b + tsing). 


Then by the Chain Rule, 


dz Ozdx  0zdy Oz Oz. 
- = —cos% + =—sina. 
x oy 


IABP) on ay By ato 


EXAMPLE 2. Find the gradient of f(x, y) = xy + y? and use it to find the derivative 


in the direction of 
ef deck ihe 
fds y=y, f(x,y) = x + 2y. 
grad f(x, y) = yi + (x + 2y)j. 


J3 J3 


re en 
2 +i 


Jul, ¥) = 


I 1 
) _ 2 = 
y+ at y) 5x + 5 


We can use Theorem | to give a geometric interpretation of the gradient 
vector. Let us assume that f(x, y) is smooth at a point (a, b), and see what happens 
to the directional derivatives f,(a,b) as the unit vector U varies. If both partial 
derivatives f,(a, b) and f,(a, b) are zero, then the gradient vector and hence all the 
directional derivatives are zero. Suppose the partial derivatives are not both zero, 
whence grad f 4 0. Then 


ju = U~+grad f = |grad f| cos 0 


where @ is the angle between U and grad /. Therefore fj, isa maximum when cos 0 = | 
and 0 = 0, a minimum when cos 6 = —1 and @ = n, and zero when cos 0 = 0 and 
0 = n/2. We have proved the following corollary. 


COROLLARY 1 


Suppose z = f(x, y) is smooth and grad f # 0 at (a,b). Then the length of 
grad f is the largest directional derivative of f, and the direction of grad f 
is the direction of the largest directional derivative of f. 


On a surface z = f(x, y), the direction of the gradient vector is called the 
direction of steepest ascent, and the direction opposite the gradient vector is called 
the direction of steepest descent (Figure 13.1.4). 


COROLLARY 2 


Suppose z = f(x, y) is smooth and 62/dy # 0 at (a,b). Then grad f is normal 
(perpendicular) to the level curve at (a,b). That is, grad f is perpendicular 
to the tangent line of the level curve (Figure 13.1.5). 
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iy 


steepest ascent 


Figure 13.1.4 


Figure 13.1.5 


PROOF By the Implicit Function Theorem, the level curve 
f(x, y) — f(a, b) =0 
has the tangent line 


Oz, i oe 
Ax at ay b) = 0. 


(a, b) is on this line. Let (xg, yo) be any other point on the line. Then 
D = (Xp — a)i + (Yo — Dj 


is a direction vector of the line, and 


6z Oz 
D- grad f = (x9 — a) + (Yo ane = 0. 


Thus grad f is perpendicular to the direction vector D. 


Water always flows down a hill in the direction of steepest descent. Thus 
on a topographic map, the course of a river must always be perpendicular to the 
level curves, as in Figure 13.1.6. 
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Figure 13.1.6 


EXAMPLE 3 A ball is placed at rest on the surface z = 2x? — 3y* at the point 
(2, 1, 5) (Figure 13.1.7), Which direction will the ball roll? 


The ball will roll in direction of steepest descent, given by —grad z. 


Oz Z 

dz =—i+ —j = 4xi — 6yj = 8i — 6j. 
grad z = —i ay xi — 6yj j 
—gradz = —8i + 6j. 


Figure 13.1.7 


The unit vector in this direction is 


—8i + 6j 8, 
= =—-—1+ 
£18? a6? 10 


Directional derivatives and gradients for functions of three variables are 
similar to the case of two variables. 
DEFINITION 


Given a real function w = f(x, y, z) and a unit vector 


U = cos ai + cos fj + cos yk 


13.1 DIRECTIONAL DERIVATIVES AND GRADIENTS 


in space, the derivative of f in the direction U and the gradient of f at (a, b,c) 
are defined as follows. 
f(a + tcosa,b + tcos B,c + tcosy) — f(a, b,c) 

t + | 


Sula, b, c) = lim 
t>0 
Ow 


ae OW + ows + —-k 
ee Be aye pe 


THEOREM 2 


Suppose w = f(x, y, Z) is smooth at (a, b,c). Then for any unit vector 
U = cos ai + cos fj + cos yk, 
the directional derivative f(a, b, c) exists and 


dw 


cos B + hi OS 
By COS y. 


Ow 
Sula, b, c) = U-grad f = By COS + az 


Corollaries 1 and 2 also hold for functions of three variables. In Corollary 2, 
grad f is normal to the tangent plane of the level surface f(x, y,z) — f(a, b,c) = 0 
at (a, b, c). 


EXAMPLE 4 Given the function 
w= zcosx + zsiny 


at the point (0, 0, 3), find the gradient vector and the derivative in the direction 


of 
1, 2 
ow Ow 
Hip ea ea 
grad w ae aye dae 


= —zsin xi + zcos yj + (cos x + sin y)k 


= —3-sin 0i + 3-cos Oj + (cos 0 + sin O)k 


= 3j+k. 
fu(0, 0, 3) = U- grad w = =+0 — 73 + el = -5. 
EXAMPLE 5 Find a unit vector normal to the surface 
z=x?4+2y?4+1 
at (1, 2, 10) shown in Figure 13.1.8. 
Let f(xy, 2) = —z + x? 4 2y? 41. 


By Corollary 2, grad f is normal to the given surface —z + x? + 2y?+1=0. 
We compute 


grad f = 2xi + 4yj — k. 
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2! (1, 2, 10) 


Figure 13.1.8 


At (1,2, 10), grad f = 2i + 8j — k. The required unit vector is found by 
dividing grad f by its length, 


2+8j)-k 248j)-k 


J? +8? +2 /69 


PROBLEMS FOR SECTION 13.1 


In Problems {-14, find the gradient vector, grad /, and the directional derivative f,. 


1 f(xy =x +97, UR 7 
2 fmy=xrty, U= aA 
3 T(x, y) a x?y3 We 31 = 4j 
4 f(xy = xy, Ua a t 4j 
Site 
5 f(x, y) = cos x sin y, ad 
ai + bj 
6 f(xy =e™*, Us ne 
a? + b? 
7 fie aa/x? + y*, v= 
2 
8 Heys fe y, U 4i 3j 
j-—2 
9 fs 9.2) = xyz, U=+ as “A : 
: Ba G2 gs 22 _i+tjt+k 
10 {m%yAJ=xr+y? 42, Uz ie 


13.2 LINE INTEGRALS 
see 2a. i-k 
11 POST tz U W 
Ps pee i-j+k 
x yp 2 J3 
13 f(% yz) = ./x? + y? +27, U =cosai + cos fj + cos yk 
14 f(x, y,2z) = Ax + By + Cz, U =cosai + cos fj + cos yk 
15 Find the derivative of z = Jxly at the point (1, 1) in the direction U = (i — 3)//10. 
16 Find the derivative of z = 1/(x + y) at the point (2, 3) in the direction U = (—i — j//2. 
17 Find the derivative of z = 2x? + xy — y? at the point (2, 1) in the direction U = ai + bj. 
18 Find the derivative of w = ./xyz at (1, 1, 1) in the direction U = (21 + j + 2K)/3. 
19 Find the derivative of w 4 — x? — y? — 2? at (1,1, 1) in the direction U = 
(i — j+ W//3. 
20 Find the direction of steepest ascent on the surface z = 2x? + 3y? at the point (1, — 1). 
21 Find the direction of steepest descent on the surface z = ./4 — x? — y? at the point 
(1, 1). 
22 Find a unit vector normal to the sphere x? + y? + z? = 9 at the point (1, 2, 2). 
23 Find a unit vector normal to the ellipsoid 4x? + y? + 42? = 3 at the point (2, 1, 3). 
O 24 Given a unit vector U = ai + bj and a function z = f(x, y) with continuous second 


partials, find a formula for the second directional derivative fyy(x, y), i.e. the derivative 
ot fy(x, y) in the direction U. 


oO 25 Given unit vectors U = u,i+ uj and V = v,i + vj, and a function z = f(x, y) with 
continuous second partials, find a formula for the mixed second directional derivative 
(fulv(%, ¥). 


13.2. LINE INTEGRALS 


There are two ways to generalize the integral to functions of two or more variables. 
One way is the line integral, which we shall study in this section. The other is the 
multiple integral, which was studied in Chapter 12. 

The line integral can be motivated by the notion of work in physics. The 
work done by a constant force vector F acting along a directed line segment from 
A to B is the inner product 


W=F-S 


where S is the vector from A to B (Figure 13.2.1). 


Figure 13.2.1 
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If the force vector 
F(x, y) = P(x, yi + O(x, y)j 


varies with x and y and acts along a curve C instead of a straight line S, the work 
turns out to be the line integral of F along the curve C (Figure 13.2.2). 


Figure 13.2.2 


The intuitive idea of the line integral is an integral 
W= [ F-dS 
c 


of infinitesimal bits of work 
dw=F-dS 


along infinitesimal pieces dS of the curve C. We now give a precise definition. 
An open rectangle is a region of the plane of the form 


ay < xX <4, by <y< by 


where the a’s and b’s are either real numbers or infinity symbols (Figure 13.2.3). 


Loo 


+—~ 


| 
I 
t 
ay ada 


An open rectangle 
Figure 13.2.3 


13.2 LINE INTEGRALS 
DEFINITION 


A smooth curve from A to B is a curve C given by parametric equations 
x = g(s) y = Als) O<s<L, 
where: A = (g(0), (0), B = (g(L), AL), 
L = length of curve, 
s = length of the curve from A to (x, y), 
dx/ds and dy/ds are continuous forO <s < L. 


We call A the initial point and B the terminal point of C. A smooth curve from 
A to B is also called a directed curve, and is drawn with arrows. 
Given s and an infinitesimal change As = ds, we let, 
Ax = g(s + As) — g(s), dx = g'(s) ds, 
Ay = h(s + As) — h(s), dy = hs) ds, 
AS = Axi + Ayj, dS = dxi + dyj. 


Thus AS is the vector from the point (x, y) to (x + Ax, y + Ay) on C, and dS is an 
infinitesimal vector tangent to C at (x, y) (Figure 13.2.4). 


y 
B 
Cc 
A 
x 
A smooth curve 
Figure 13.2.4 
DEFINITION 
Let F(x, y) = P(x, yi + Q(x, y)j 


be a continuous vector valued function on an open rectangle D and let C bea 
smooth curve in D. The line integral of ¥ along C, 


[Fas | Pdx + Qdy, 
c _ vc 


is defined as the definite integral 
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Notice that the inner product of F and dS is 
F-dS = (Pi + Qj)-(dxi + dyj) = Pdx + Qdy. 


This is why we use both notations |. F+ dS and {- P dx + Q dy for the line integral. 


DEFINITION 


The work done by a continuous force vector F(x, y) along a smooth curve C 
is given by the line integral 


w= | F-dS. 
c 


JUSTIFICATION We can justify this definition by using the Infinite Sum Theorem 
from Chapter 6. Let W(u, v) be the work done along C from s = u tos =v 
(Figure 13.2.5). Then W(u, v) has the Addition Property, because the work 
done from u to v plus the work done from v to wis the work done from u to w. 
On an infinitesimal piece of C from s to s + As, the work done is 


AW = F(x, y)- AS & F(x, y)« dS (compared to As). 


dx 
But F(x, y)+dS = Pdx+ Qdy= tes + of ds. 


Figure 13.2.5 


By the Infinite Sum Theorem, 
Bl dx dy 
Ws P ds = | F-dS. 
i) | ds +02 ‘ if 


The next theorem is useful for evaluating line integrals. It shows that any 
other parameter t can be used in place of the length s of the curve. Figure 13.2.6 
illustrates the four parts of this theorem. 


13.2 LINE INTEGRALS 


A(Xo, Yo) 


» y B(xo, yi) 
F 
A(X, Yo) y B(x1, Yo) 
F 
x 


(i) Horizontal (ii) Vertical 


(iii) Parametric curve 


y 


J, F-ds — [, F-ds 


(iv) Reversing the curve direction 
Figure 13.2.6 


THEOREM 


Let | F - dS be a line integral. 


(i) If C is a horizontal directed line segment xy S x S X1,yY = Yo, then 


| F-dS = P(x, Yo) dx. 
Cc xo 


(ii) If C is a vertical directed line segment x = X9,Vo S y S yy, then 


yt 
[ Peas = | OG0,.n0ds. 
c yo 
(iii) If C is traced by a parametric curve 
x=g(), y=hit), e<tse 


where dx/dt and dy/dt are continuous, then 
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[F dS = {. (PG+0 a a 


(iv) Reversing the curve direction changes the sign of a line integral. That is, 
if C, is the curve C with its direction reversed, then 


{ F-dS = -| F-dS. 
C1 c 


Xy yi 
Remark The integrals | P(x, Yo) dx, | OAxXo, y) dy 


are sometimes called partial integrals. 


PROOF (i) and (ii) are special cases of (iii). (ili) is proved by a change of variables, 
d 
{¥ dS = ‘ (Po: +0 i a 
dx dy ds 
(Pa + ala 
= i. [Pa + of at. 


(iv) is true because reversing the limits changes the sign of an ordinary 
integral. 


EXAMPLE 1 Find the line integral of 
F(x, y) = sin x cos yi + ej 
along 
(a) The horizontal line C,:0 < x < 2, y = n/3 (Figure 13.2.7(a)). 
(b) The vertical line C,:0 < y < 1,x = 2 (Figure 13.2.7(b)). 


Figure 13.2.7 (a) (b) 


We use partial integrals. 


(a) { Feds = [ sin x cos ~ dx 
Cy 0 3 


= Isin dx = —Leos.| Sohal eat 


0) ie) 
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1 1 
(b) | F-dS= { e” dy = vad = He? — 1). 
0 


C2 i) 


Given two points A and B, there are infinitely many different smooth curves 
C from A to B. In general the value of a line integral will be different for different 
curves from A to B. 


EXAMPLE 2 Let the force vector F be F = — yi + xj. 
F is perpendicular to the position vector xi + yj but has the same length 
as xi + yj. Find the work done by F along the following curves, shown in 
Figure 13.2.8, from (0, 0) to (1, 1): 
(a) C,:Theliney=x,0Sx <1. 
(b) C,: The parabola y = x*,0< x <1. 
(c) C3: The curve y = V/x,0<x< 1. 


y 
dl. 1) 
Cy 
Oo X 
(a) 
Figure 13.2.8 


(a) Putx=t,y=t. 


w= | F-dS= —ydx+xdy 
Cy 


cy 


= [ + thdt = 0. 


The work is zero because the force F is perpendicular to dS along C,. 
(b) Putx=ty=?. 


W, -| —ydx+xdy 
C2 


1 1 
-{ (-2 + 0-29 dr = [ t? dt =+. 
0 


i) 


(c) Putx=f,y=t. 


1 
Ws = { —yde + xdy = [ (—t-3t? + A)dt 
: C3 (0) 


1 
=| 207 dt = —4. 


¢) 
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A piecewise smooth curve is a curve C that can be broken into finitely many 
smooth pieces C;, C2,..-, C, where the terminal point of one piece is the initial 
point of the next (Figure 13.2.9). For example, a curve formed by two or more sides 
of a rectangle or a polygon is piecewise smooth. The line integral of F(x, y) over a 
piecewise smooth curve C is defined as the sum 


[ Feds = F-dS + Peds +4 [ F-dS. 
6 Cy 


C2 n 


y B 
Cy 


Figure 13.2.9 A Piecewise Smooth Curve from A to B 


EXAMPLE 3 Find the line integral 


i xydx + x*ydy 
c 


where C is the rectangular curve from (2, 5) to (4, 5) to (4, 6). 


We see in Figure 13.2.10 that C is a piecewise smooth curve made up of a 
horizontal piece 


C,:2<x <4, ea) 
and a vertical piece 
Cy:x = 4, S<yS6. 


The line integral is the sum of two partial integrals, 


Figure 13.2.10 
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4 6 
[ Beas = F-as + [ Feds = [ x-Sdx + [ 4? .y dy 
c Ci C2 2 5 


4 6 
= 5-43? | + 16-4y°| = 304+ 88 = 118. 
2 5 


A simple closed curve is a piecewise smooth curve whose initial and terminal 
points are equal and that does not cross or retrace its path. Examples of simple 
closed curves are the perimeters of a circle, a triangle, and a rectangle. The value 
of a line integral around a simple closed curve C depends on whether the length s 
is measured clockwise or counterclockwise, but does not depend on the initial point 
(Figure 13.2.11). The clockwise and counterclockwise line integrals of F around a 
simple closed curve C are denoted by 


Feds, f Fas. 
Cc Cc 


Q F-ds @ F-ds 
C F F 
p F-dS Q F-d8 


Figure 13.2.11 Integrals around Simple Closed Curves 


THEOREM 2 


If C is a simple closed curve, then 


> F-dS = — > Fas 
c Cc 


and the values do not depend on the initial point of C. 


PROOF The equation in Theorem 2 holds because reversing the direction of the 
curve changes the sign of the line integral. Suppose C has the initial point A, 
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and its direction is clockwise. Let A, be any other point on C, and let C, 
and C, be as in Figure 13.2.12. 


With the initial point A, 


f Fd = FedS + F-dS. 
Cc Cy 


C2 


With the initial point A,, 


Peas = [ P-dS + [ F-dS. 
Cc C2 Ci 


These are equal as required. 


Figure 13.2.12 


EXAMPLE 4 Find the line integral 


> —ydx +xdy 
Cc 


where C is the circle x? + y* = 4, shown in Figure 13.2.13. 


(2,0) * 


Figure 13.2.13 


We may start at any point of C. Take (2, 0) as the initial point. Then C has 
the parametric equations 


x = 2cos 0, y = 2sin 0, 0<0< 2n. 


As @ goes from 0 to 2z, (x,y) goes around C once counterclockwise as 
required. 
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an dx dy 
$ —ydx+xdy =| [— 955 + = dé 


2n 
a (—2 sin 6(—2 sin #) + 2 cos 6(2 cos €)) dé 


¢) 


2n 2n 
= 4sin?0 + 4cos? od0 = | 4d9 = 8n. 


0) 0 


Line integrals in space are developed in a similar way. Instead of an open 
rectangle we work in an open rectangular solid. A smooth curve C in space has three 
parametric equations with continuous derivatives, 


x = 2(s), y = h{s), z = ([(s), O<s<L. 
Given a continuous vector valued function 
F(x, y, z) = P(x, y, z)i + O(x, y, z)j + R(x, y, z)k 


and a smooth curve C in space, we define the line integral of F along C, in symbols, 


[ F-as= | Pdx + Qdy + Rdz, 
c c 


Lf dx dy dz 
as [ te + QT + nea. 


EXAMPLE 5 Find the line integral 
Zz 
{ (x + y)dx + —dy + xydz 
Cc x 
along the spiral C given by 
x = cost, = sint, 7 = 2t, O<t< 
The line integral is 


mi? 2t 
i) (cos t + sin t) d(cos t) + —— d(sin t) + (cos t sin t) d(2t) 
3 cos t 


{2 
= [ (—costsint — sin*t + 2t + 2costsin t) dt 
0 

n{2 


2 


Le 2 1 ae 2 2 
= —~sin*t — |~f — =sintcost}) + t° + sin‘t 
2 <9 é 


dati Poh e gag st + t? ee 
=-simMn = =~ $n? CO. 
y) pak) ‘ 


rel 
4. 
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PROBLEMS FOR SECTION 13.2 


Evaluate the following line integrals. 


1 


19 


20 


{ xe dx + x*ydy, C:0<x<2,y=3 
Cc 
{ x dx + x*ydy, C:0<y<4,x=4 
c 


| xe dx + x*ydy, C:x=3ty=?,0<t<1 
c 


[_serdx + sy dy, C:x=evp=e,-lLst<l 
c 


{ (00s xi + sin yf) -8, C:x=ty=t05r<sl 
Cc 


{ & + J | -dS, Cis the rectangular curve from (1, 1) to (3, 1) to (3, 6). 
c\xXy N+P) 


[lat Las, C:x =2,y=51<1<4 
fa xy xX } 


{ fa J | as, Cix=ty=P,l<is4 
c\xy x+y 

$ ydx —xdy anddp y dx —xdy, C:ix*+y?=1 

Cc 

$ xeydx + xy? dy, Cire +y?=4 

Cc 

¢ (x + y)dx — 3xydy, C:x? +)? =4 

Cc 

$ te cos yi + e*sin yj)+dS, C is the square with vertices (0, 0), (1,0), (1, 1), (0, 1). 


f (/xyi + x?y?j)-dS,_ C is the triangle with vertices (0, 0), (1, 1), (1, 0). 
Cc 


{ yzdx + xzdy+xydz, C:x=ty=e.z=P,0<t< 1, 
Cc 


i yedx + xzdy + xydz. C:x =cost.y =sint,z = tant,0<¢< n/4, 
Cc 


[ (xi + yj + 2k)+ dS, C is the rectangular curve from (0, 0,0) to (1, 0,0) to (1, 1,0) 
Cc 
to (1, 1, 1). 


Find the work done by the force F = (xi + yj)/(x? + *) acting along a straight line 
from (1,1) to (2,5). 

Find the work done by the force F = (iy + 1)) — (jAx + 1)) acting along the parabola 
x=hy= rocr<i. 

Find the work done by the force F = x7i + y?j + z°k acting along a straight line from 
(0, 0, 0) to (3, 6, 10). 

Find the work done by the force F = yi + zj + xk along the curve x = sf ts yo W/t, 
z=t1<r1<4. 


13.3. INDEPENDENCE OF PATH 


13.3 INDEPENDENCE OF PATH 


For functions of one variable, the Fundamental Theorem of Calculus shows that 
the integral is the opposite of the derivative. In this section we shall see that the line 
integral is the opposite of the gradient. 

By a vector field we mean a vector valued function 


F(x, y) = P(x, yji + O(x, y)j 


where P and Q are smooth functions on an open rectangle D. 

For example, if f(x, y) has continuous second partials on D then its gradient 
grad f is a vector field. 

Many vector fields are found in physics. Examples are gravitational force 
fields and magnetic force fields, in which a force vector F(x, y) is associated with each 
point (x, y). Another example is the flow velocity V(x, y) of a fluid. A vector field in 
economics is the demand vector 


D(x, y) — D(x, y)i + D,(x, y)j, 


where D,(x, y) is the demand for commodity one and D,(x, y) is the demand for 
commodity two at the prices x for commodity one and y for commodity two. All of 
the examples above have analogues for three variables and three dimensions (and 
the demand vector for n commodities has n variables and n dimensions). 


DEFINITION 


I (x, y) is a potential function of the vector field Pi + Qj if the gradient of f is 
Pi+ Qj. 


Not every vector field has a potential function. Theorem 1 below shows 
which vector fields have potential functions, and Theorem 2 tells how to find a 
potential function when there is one. 

Using the equality of mixed partials, we see that if the vector field Pi + Qj 
has a potential function, then 0P/dy = 0Q/0x. If fis a potential function of Pi + Qj, 
we have 

0 ) 
grad f= Ti+ Ti = Pi + Ob 
OP Of =—rf SD 
dy dyéx Odxdy ex’ 


EXAMPLE 1_ The vector field — yi + xj has no potential function, because 


OP a@&-y) 1 6Q Ox | i 
dy ey °° ox ax 
THEOREM 1 
: ’ ; fates _ OP dO 
A vector field Pi + Qj has a potential function if and only if oy Tae 
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We have already proved one direction. We postpone the proof of the other 
direction until later. 

By definition, grad f = Pi + Qj if and only if df = P dx + Q dy. In general, 
an expression P dx + Q dy is called a differential form. A differential form is called 
an exact differential if it is equal to the total differential df of some function f(x, y). 

Using this terminology, Theorem | states that: Pdx + Qdy is an exact 
differential if and only if @P/éy = 0Q/dx. 


EXAMPLE 2 Test for existence of a potential function: 
x*yi + sin x cos yj. 
éP A(x’y) 

éy say 


-2 


x, 


6Q- Asin x cos y) 
—= = = COS X COS y. 
ex ex 


There is no potential function. 


EXAMPLE 3 Test for existence of a potential function: 
3x7 y7i + (y? + 2x3 y)j. 
OP _ A3x"y") _ 62 


a x ys 

oy dy 

AQ Ay? + 2x3y) | 6x2 
ax ex Er 


There is a potential function. 


THEOREM 2 (Path Independence Theorem) 


Let Pi + Qj be a vector field such that @P/dy = 0Q/0x and let A and B be 
two points of D. 


(i) Let f be a potential function for Pi + Qj. For any piecewise smooth 
curve C from A to B, 


{ P dx + Qdy = f(B) — f(A). 
Cc 


Since the line integral in this case depends only on the points A and B and 
not on the curve C (Figure 13.3.1), we write 


B 
i Pdx+ Qdy= i Pdx + Qdy. 
A c 
(ii) g is @ potential function for Pi + Qj if and only if g has the form 


(xy) 
g(x, y) = i Pdx+Qdy+kK 


A 


for some constant K. 
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y 


[UP dx+Qdy 


Figure 13.3.1 Independence of path 


Theorem 2 is important in physics. A vector field of forces which has a 
potential function is called a conservative force field. The negative of a potential 
function for a conservative force field is called a potential energy function. Gravity, 
static electricity, and magnetism are conservative force fields. Part (i) of the theorem 
shows that the work done by a conservative force field along a curve depends only 
on the initial and terminal points of the curve and is equal to the decrease in potential 
energy. 

Mathematically, Theorem 2 is like the Fundamental Theorem of Calculus. 
It shows that the line integral of grad f along any curve from A to B is equal to the 
change in the value of f from A to B. When A = B, we have an interesting con- 
sequence: 


If f(x, y) has continuous second partials then the line integral of the gradient 
of f around a simple closed curve is zero, 


$ grad f-dS = 0. 
Cc 


Using part (ii), we can find a potential function f(x, y) for a vector field 
Pi + Qj in three steps. 


oP 0d 
== D. 
When to Use dy Ag on 


Step 7 Choose an initial point A(a, b) in D. 


Step 2 Choose and sketch a piecewise smooth curve C from A to an arbitrary 
point X(xXo, Yo). 


Step 3 Compute f(X9, Yo) by evaluating the line integral 
flo. yo) = | Pax + Ody. 
c 


We postpone the proof of Theorem 2 to the end of this section. 


EXAMPLE 3 (Continued) Find a potential function for the vector field 


3x7y*i + (y? + 2x3 y)j. 
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OP 
We have already shown that — = og 

oy @x 
Step 7 Pick (0, 0) for the initial point. 


Step 2. Let C be the rectangular curve from (0,0) to (0,9) to (x9, yg), Shown in 
Figure 13.3.2. 


(0, 39) (Xo, Yo) 


(0, 0) 


Figure 13.3.2 


Step 3 A potential function is 


Ihe Lec= } 3x2y? dx + (y? + 2x4y) dy 
C 


yo xo 
= | (y? + 2+03y)dy + { 3x? dx 
0 ) 
= 493 + 2898. 
f(x,y) = 4y? + xy’. 
As a check we may compute grad f. 
6 a) 
grad f = a + S = 3x7y?i + (y? + 2x3 y)j. 
Ox oy 
We can get the same answer by choosing another curve in Step 2. 


FIRST ALTERNATE SOLUTION 


Step 2 Let C, be the rectangular curve from (0, 0) to (x9, 0) to (Xo, Yo), shown in 
Figure 13.3.3. 


(xo, Yo) 


(Xo, 0) 
Figure 13.3.3 
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Step 3 fl%o Yo) = { 3x2y? dx + (y? + 2x3y) dy 
Cy 


xO yo 
=| 3x? *0? dx + | (y? + 2x3y) dy 
Jo 0 
= 0 + Gyo + xoyo)s 
L(x, y) = Sy? + x? y?. 
SECOND ALTERNATE SOLUTION 


Step 2 Let C, be the straight line from (0,0) to (x9, yo), shown in Figure 13.3.4. 
It has parametric equations 


xX = [Xq, y=tyo, O<t<l. 


(Xo, Yo) 
Co 


(0, 0) x 
Figure 13.3.4 


Step 3 f(xX9,Vo) = | 3x?y? dx + (y? + 2x4y) dy 


C2 


1 
= { [3(tx9)*(t¥o)2}xo + [ltvo)? + 2Utxo)?(t¥o)]%o dt 
ré) 


1 
= { 3xdyt* + t7y38 + 2t*xhy6 dt 
C8) 


Fe ey antes tor PS ee ee 
= $X0Vo + 3Yo + SXoVO = XoVo + 30, 


f(x, y) = xy? + Sy? 


EXAMPLE 4 An object at the origin (0,0) has a gravity force field with magnitude 
proportional to 1/(x? + y?) and the direction of —xi — yj. Show that this 
force field is conservative and find a potential function. 


The force vector is 


—xi- yj) k 
Ft 9) =| | 2 2 
/x? 4 y? x +y 


= —kx(x? + y?)79?i — kyl? + y?) ?j, 


for some constant k. F(x, y) is undefined at (0, 0) but is a vector field on the 
open rectangle 0 < x. 
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oP = ~kst){ - 3) + y2)7 8/2 = 3kxy(x? + y2)7 9, 
oy 2 
20 _ key) —3] 02 + y2)-9? = Bhp? + 87. 
ox 2 


Therefore F is conservative. 
Step 7 Take the initial point (1, 0). 


Step 2 Let C be the rectangular curve from (1,0) to (1, yo) to (Xo, Yo), Shown in 
Figure 13.3.5. 


(1, Yo) 
(xo, ¥o) 


: (1, 0) a 
Figure 13.3.5 


yo Xo 
Step 3 f (Xo, Yo) = { —ky(l + y?)-3 dp + [ = kx(x? +b yoy 7? dx 
0 


1 


yo xo 


= kl + ft 
0 l 


= k(1 + yo)? — k + k(x2 + y2)7 2 — kL + yy?) 1? 
= k(x2 + y2)"1/? + constant, 


+ K(x? + aie | 


T(x y= a + constant. 


Jx? + y? 


Any choice of the.constant will give a potential function. The same method 
works on the open rectangle x < 0. 


An exact differential equation is an equation of the form 
P(x, y)dx + Q(x, y)dy = 0, 


where 6P/dy = 0Q/dx. Exact differential equations can be solved using Theorem 2. 


EXAMPLE 5 Solve the differential equation 
(x? + sin y)dx + (x + 1)cos ydy = 0. 
First we test for exactness. 


A(x? + sin y) a((x + 1) cos y) 
Oy, oe ee ae, a OS 


Next we find a function with the given total differential. That is, we find a 
potential function for the vector field 


(x? + sin y)i + (x + 1)cos yj. 
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Step 7 Take (0, 0) for the initial point. 


Step 2. Let C be the rectangular curve from (0,0) to (0, yo) to (x9, Yo), Shown in 
Figure 13.3.6, 


(0, Yo) (Xo, Yo) 


(0, 0) x 
Figure 13.3.6 


yo xo 
Step 3 f (Xo; Yo) = | (0+ 1)cos ydy + [ (x? + sin yo) dx 
1) 0 
sin Yo + 3X + Xo SiN Yo, 
f(x,y) = siny + 4x7 4 xsiny. 
Step 4 f(x,y) isa constant k because df = 0. The general solution is 


Ul 


4x3 4+ xsiny + siny = k, 
sin aoa 
cea Tas 


ais 


: 3 
= arcsin k constant. 
y x+1’ 


We conclude this section with the proofs of Theorems 1 and 2. The proof 
of Theorem 1 uses a lemma about derivatives of partial integrals. 


LEMMA 


Suppose P(x, y) is smooth on an open rectangle containing the point (a, b). Then 


<{ Plt, y) dt = P(x, y), 
Ox Ja 


a px * AP 
=| P(t, y)dt = | ee. y) dt. 
OY Ja a OY 


PROOF The first formula follows at once from the Fundamental Theorem of 
Calculus. For the second formula, let Ay be a nonzero infinitesimal and let 


z= ) P(t, y) dt. 
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As y changes to y + Ay, we have 


x 


[ " P(t,y + Aydt — i; P(t, y) dt 
Az dg 


a 


Ay Ay 
* P — P(t * Op 
| (t, y + Ay) POLY) = | ands 
i Ay a OY 
qotnewenasapen: ~e26 EG a 
aking Stan ar par ; ay Ai +y¥ i 


PROOF OF THEOREM 71 We must find a potential function for Pi + Qj. 
Assume 0P/d) = 6Q/dx. Pick a point (a, b) in D, and let f(x, }9) be the line 
integral of Pi + Qjon the rectangular curve C from (a, b) to (a, yp) to (x9, Yo) 
(Figure 13.3.7). Thus 
flo. ¥) = | Pax + Qdy 
Cc 


yo xo 
a { O(a, Yo) dy + { P(x, Yo) dx. 
b a 


(a, Ya) (Xo, ¥o) 


Figure 13.3.7 


By the Lemma, 


gf 7 sll Oa) aye iQ P(x, y) ix = P(x, y). 
= an O(a yay + [Pls ya 
= O(a.) + [ Found dx 
= ota + [Sos nar 
= O(a, ») + (Ol, ») - Ola, y)] = Q0e,). 
Thus Farey, F= o6,y, 


and df = Pdx + Qady. 
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PROOF OF THEOREM 2 


(i) Let C have the parametric equations 


(ii) 


Theorems 1 and 2 also hold for three variables. For three variables a vector 


x = a(t), y = h(0), C, St<ec. 
Then A = (g(c,), h(c,)) and B = (g(c,), h(c)). By the Chain Rule, 


dz _ 0zdx | dz dy 
dt dxdt dydt 


= Pg'(t) + Oh'(t). 


Then | P(x, y)dx + Q(x, y)dy = [- Pe'(t) + Oh'(t) dt 
Cc 


cy 


= f(g(c2), A(c2)) — F(gler), hes) 
= f(B) — f(A). 


Asimilar computation works for piecewise smooth curves. This proves (i). 


Define f(x, y) by 
x 
f(X) = { P dx + Qdy, 
A 


where A = (a,b), X = (x, y). Let C be the rectangular curve from (a, b) 


to (a, y) to (x, y). Then 


Stsy= | Pax + QO dy. 


We already showed in the proof of Theorem 1 that this function 
f(x, y) is a potential function for Pi + Qj. To complete the proof we 


note that the following are equivalent. 


grad g = Pi + Qj, 


ee og _ of 

6x Ox dy ay’ 

OCA) 2G eg OCI) Le 
Ox oy 


g — f depends only on y and only on x, 
a(x, y) v? F(x, y) = constant, 


xX 
g(x, y) = | Pdx + @dy + constant. 
A 


field has the form 


F(x, y, z) = P(x, y, z)i + Q(x, y, z)j + R(X, y, zk. 


Theorem 1 for three variables reads as follows. 
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THEOREM 1 (Three Variables) 


A vector field Pi + Qj + Rk has a potential function if and only if 


6P_@Q OP OR 0Q_ OR 
dy ax’ dz Ox’ éz ay” 


Theorem 2 is modified in the same way. 


PROBLEMS FOR SECTION 13.3 


Test the following vector fields for existence of a potential function and find the potential function 
when there is one. 


Corn aun & Ww NY — 


=—=— — 
nH = & 


— 
Lv) 


— ee 
n vw BS 


v= = 
co Oo @ ~I 


21 


(2x + y?)i + (x? + 2y)j 

im yA 

yi + 2xj 

xei t+ yery 

Jx? + yi+j) (x >0,y > 0) 
y cos xi + ysin xj 

y cos xi + sin xj 

eri + j) 

—2i + 6j 

ple yitx/x?t yj (x>0y>0) 
xpi t+ xy4j 

447) &>Oy>0) 

xy 

(3x + Sy)i + (Sx — 2y)j 


y? 
{+ 2ylnxj (x > 0) 
x 


sinh x cosh yi + cosh x sinh yj 


JSylxi + /x/yj_ (x > 0, > 0) 

Show that every vector field of the form P(x)i + Q(y)j has a potential function. 

Show that every vector field of the form f(x + y)(i + j) has a potential function. 

Show that every vector field of the form f(x? + y?)(xi + yj) has a potential function. 
Show that every vector field of the form f(xy)(yi + xj) has a potential function. 


Show that the sum of two conservative force fields is conservative. 


In Problems 22-31 solve the given exact differential equation. 


22 
23 
24 
25 


e*dx + sinydy =0 
(3x + 4y)dx + (4x — 2y)dy = 0 
(x? + 2xy + y*)dx + (x? + 2xy + 8) dy =0 


(Jx + /y) dx + (x/2/y) dy =0 
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26 2x sin ydx + (y + x? cos y)dy =0 
y 2 = 
27 Ris {eto + arctan x)dy = 0 
28 (ax + by) dx + (bx + cy)dy = 0 
29 sin x sin ydx — cosxcosydy =0 


arcsin y ed: In x 
Ji-y 
31 (x+ ./x+y)dx+(y+/x+ y)dy=0 


dy=0 


32 Find a function Q(x, y) such that Jxy? dx + Q(x, y) dy is an exact differential. 

33 Find a function P(x, y) such that P(x, y) dx + sin? x cos y dy is an exact differential. 

34 The gravity force field of a point mass in three dimensions has magnitude proportional 
to i/(x? + y* + 2) and the direction of —xi — yj — zk. Show that the force field is 
conservative. 
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Green’s Theorem gives a relationship between double integrals and line integrals. 
It is a two-dimensional analogue of the Fundamental Theorem of Calculus, 


b 
F(b) — F(a) = { F'(x) dx, 


a 


and shows that the line integral of F(x, y) around the boundary of a plane region D 
is equal to a certain double integral over D. 
Let D be a plane region 


4;5xXS4), bi(x) < y S bo). 


The directed curve which goes around the boundary of D in the counterclockwise 
direction is denoted by 0D and is called the boundary of D (Figure 13.4.1). 


y 


0D 


Figure 13.4.1 The boundary of D 


If b, (x) and b,(x) have continuous derivatives, dD will be a piecewise smooth 
curve and thus a simple closed curve (see Section 13.1). 


GREEN‘’S THEOREM 


Let P(x, y) and Q(x, y) be smooth functions on a region D with a piecewise 
smooth boundary. Then : 
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Pd dy = ea 
f., x + Qdy IVs ayo 


¢. —Qdx+Pdy= ie + Baa. 


(See Figure 13.4.2.) 


Figure 13.4.2 a 


The second formula follows at once from the first formula by replacing 
P by —Q and Q by P. We shall prove the theorem only in the simplest case, where 
D is a rectangle. 


PROOF FOR D A RECTANGLE D is shown in Figure 13.4.3. 


¥ 
by t--~-- 
betes I 1 
i} I 
a t ie 
a, ay 


Figure 13.4.3 


The line integral around @D is a sum of four partial integrals 


ar Da 
f Pax + Qdy=[ Pesbydx + [ Qlar.sar 
cD a1 m 


+ 


ay by 
+ | P(x, by) dx + { O(a,, y)dy 
a2 bs 


be a> 
=[ Otar. 9) - Oa, yay — f° Pex.ba) ~ Pox bi) dx 
hy ay 


By the Fundamental Theorem of Calculus, 


ta a 

QO(a2, y) — Oa,,y) = [ ae 1x, 
he 6 

P(x, by) — P(x, by) = aes 
by oy 
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Therefore 


b> To 22 bo oP 
$ Pdx + Qdy = [ i ax dy ~ [ i ao 
6D bi Yai ay 1 
az pba oP 
-[ fe cr axdy = [[ 2 ae 
a, Jb, OX OY E Cx oy 


We may apply Green’s theorem to evaluate a line integral by double 
integration, or to evaluate a double integral by line integration. 


EXAMPLE 1 Compute the line integral 


$ pee ay 
6D 


by Green’s Theorem, where D is the rectangle shown in Figure 13.4.4, 


O<x<2, O<sysl. 


Figure 13.4.4 


By Green’s Theorem, 


Ax + y)  A(x?y) 
5 xpd dy = dA 
$y x + (x + y)dy {J ax 


dy 


= {fe —x*)dA= { f 1 — x? dy dx 
2 : 


2 
={ 1 —x?dx = —2. 


(¢) 


As a check, we also compute the line integral directly. 


2 1 
29 dx + + ydy = [ 20x + | 2+ ydy 


10) 0 


0 0 
+f x? +1dx +| 0+ ydy 
2 1 


58 4 2 
=0+3-3-7= -3. 


EXAMPLE 2 Evaluate by Green’s Theorem the line integral 


y 
——dx +2 
> ‘haal xy dy 
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where D is the region bounded by the curve y = x* and the line y = x, 
shown in Figure 13.4.5. 


y 
aD 
x 
Figure 13.4.5 
D is the region Oza x Ya ype yy, 
, (2 1 
6 y “dx + xy dy ars xy) ayix +) a, 
ep éy 
NL Opes 
{fa rae 
D 
1 px 
=| 2y — ——~ dy dx 
0 vx 
i 2 
=| es ogee eee a 
0 x+1 x41 
=2In2-% 


As a corollary to Green’s Theorem we get a formula for the area of D. 


COROLLARY 


If D has a piecewise smooth boundary, then the area of D is 


A= > xdy = > — ydx. 
6D 6D 
PROOF By Green’s Theorem, 


$ r= [15 - 544 = [faana 
fy [f aa [Jana 


EXAMPLE 3 Use Green’s Theorem to find the area of the ellipse shown in Figure 
13.4.6, 


y 


pol 


x. 
fe 


Q 
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0D 


Figure 13.4.6 


The boundary of the ellipse is the parametric curve 
y = bsint, O<t<2z. 


x =acost, 


By the corollary, 
2n d 2n 
x at = i (a cos t)(b cos t) dt 


A=4 xdy = | 
aD 0 


2n 
= ab { cos” t dt = nab. 
oO 


Green’s theorem has a vector form which is convenient for physical applica- 
tions. We define two new functions obtained from a vector field, the curl and the 


divergence. 
DEFINITION 

Given a vector field F(x, y) = P(x, y)i + Q(x, y)j in the plane. 

The curl of F is curl F = 0 < 

ox dy 

oP 0d 

divF = — + — 

IV ax + By 


The divergence of F is 


On the boundary OD, the differential forms P dx + Q dy and —Q dx + P dy 


may be written in the vector form 
Pdx + Qdy =F-Tds, 
—Qdx + Pdy =F-Nds, 


T = unit tangent vector to dD, 


where 
T ds = dxi + dyj, 
and N = unit outward normal vector to dD, 
N ds = dyi — dxj. 


T and N are shown in Figure 13.4.7. 


820 13° VECTOR CALCULUS 


d T N 


aD 


Figure 13.4.7 


Substituting the vector notation into the original form of Green’s Theorem, 
we get the following. 
GREEN’S THEOREM (Vector Form) 


Given a vector field F(x, y) = Pi + Qj on a region D with a piecewise smooth 


boundary, 
F-Tds = iA curl F dA, 
éD 


D 


f F.Nds = ia div F dA. 
€D 
D 


The physical meaning of Green’s theorem can be explained in terms of the 
flow of a fluid (a liquid or gas). Let the vector field F(x, y) represent the rate and 
direction of fluid flow at a point (x, y) in the plane. Consider a plane region D and 
element of area AD containing (x, y) (Figure 13.4.8). 


Figure 13.4.8 


We first explain the formula 


7 F-Tds = [| cunkaa, 
“éD D 


The line integral $ F + T ds of the flow component in the direction tangent to the 
€D 


boundary is called the circulation of F around 6D. Green’s Theorem states that the 
circulation of F around the boundary of D equals the integral of the curl of F over D 
(Figure 13.4.9). 
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Figure 13.4.9 
When we apply Green’s theorem to an element of area AD we get 


eee «T ds ~ curlF AA (compared to AA). 
Thus the curl of F at (x, y) is equal to the circulation per unit area at (x, y). 

If curl F is identically zero, the fluid flow F is cailed irrotational. By the 
Exactness Criterion, F is irrotational if and only if P dx + Q dy is an exact differential. 
The circulation of an irrotational field around any OD is zero. 

Next we explain the formula 


F-Nds = [[oiveaa. 
éD 


D 


The line integral F - N ds of the flow component in the direction of the outward 
6D 

normal vector is called the flux across 0D. The flux is the net rate at which fluid is 

flowing from inside D across the boundary and is therefore equal to the rate of 

decrease of the mass inside D. Green’s Theorem states that the flux of F across the 

boundary of D equals the integral of the divergence of F over D (Figure 13.4.10). 


Figure 13.4.10 


When we apply this to AD we get 
F-Nds = divF AA (compared to AA). 


é@AD 
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Therefore the divergence of F at (x, y) is the net rate of flow of fluid away from (x, y), 
and is equal to the rate of decrease in density at (x, y). Positive divergence means that 
the density is decreasing, and negative divergence means that the density is increasing. 

If div F is identically zero, the fluid flow is called solenoidal, or incompressible. 
By the Exactness Criterion, F is incompressible if and only if —Q dx + P dy is an 
exact differential. The flux of an incompressible field across any OD is zero. 


EXAMPLE 4 A fluid is rotating about the origin with angular velocity w radians 
per second. Find the curl and divergence of the velocity field F(x, y). 


As we can see from Figure 13.4.11, the velocity at a point (x, y) is 


F(x, y) = @(—yi + xj) = —ayi + ox). 


Figure 13.4.11 


Awx) (ay) 


= 20, 


Ox oy 
—wy a 
ol-—w y) ie (ax) 


——— = 0. 


Ox oy 


Then curl F = 


divF = 


Thus a purely rotating fluid is incompressible and its curl at every point is 
equal to twice the angular velocity. 


EXAMPLE 5 A fluid is flowing directly away from the origin at a rate equal to a 
constant b times the distance from the origin (Figure 13.4.12). Find the 
curl and divergence of the flow field. 


la 
a 


Figure 13.4.12 
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We have F(x, y) = b aD ——,/x? + y? = bxi + byj. 


a 


(by)  a&bx) 


url F = Ax = ay = 0. 
dig ae) SOM) = 5p. 
Ox oy 


The fluid flow field is irrotational and the divergence at every point is 2b. 


PROBLEMS FOR SECTION 13.4 
In Problems 1-12, find the line integral by Green’s Theorem. 


1 2ydx+3xdy, D:0Sx<10sy<1 
aD 


2 $ xydx +xydy, D:O0<x<10<y<I 
aD 


3 > e**3¥ dx 4 e dy, D:-2<x<2,-1<y<1 
aD 

4 f, yeosxdx + ysin xdy, D:0<x<xn/1<y<2 

5 xiv dx + xy? dy, D:0<x<10<y<x 

6  eydx + J+ vay Dil<x<2,2xsy<4 

7 BO et Aaa ay, Dil<x<21<y<x? 

8 sin ydx + sin x dy, D:0Sx<n/2,x< y<n/2 

9 f, xin yar, Dil<x<2et<y<e” 

10 p Vit ed, D:-15x<1,x?<y<1 

11 $ xydx — xy? dy, D:x?+y? <1 Hint: Use polar coordinates. 

12 yi dx + 2x3 dy, D:x?+y? <4 


oD 


In Problems 13-18, find (a) curl F, we F. Tds, (c) divF, ¢ F-Nds. 
aD aD 


13 F(x, y) = xyi -— xyj, D:O<x<510<sy<1 

14 F(x, y) = ax*i + by?j, D:0<x<1,0Sy<1 

15 F(x, y) = ay?i + bx?j, D:0<x<1,0<y<x 

16 F(x, y) = sin x cos yi + cosxsinyj, D:0<x<27/2,0<y<x 
17 F(x, y) = yi — xj, Dix? + y? <1 


18 F(x, y) = xi+ yj, Dix?+y? <1 
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19 
20 


21 


O 22 


13.5 


23 


25 
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Use Green’s Theorem to find the area inside the curve r = a + cos @, (a = 1). 


Use Green’s Theorem to find the area inside the ellipse x?/a? + y?/b? = 1 and above 

the line y= c(0<c < b). 

Show that if D has a piecewise smooth boundary, the area of Dis A = 4 —ydx + xdy. 
oD 

Show that for any continuous function f(t) and constants a, b, ¢, 


$ af (x? + y*)dx + bf(x? + y?)dy = 0 
aD 


where D is the circle x? + y? < c?. 


Find the value of the line integral 
) (a,x + byy) dx + (a,x + bry) dy 
eD 


where D is a region with area A. 
Show that any vector field of the form 

F(x, y) = xf (x? + yi + Wf? + yj 
is irrotational. 
Show that any vector field of the form 

F(x, y) = yf? + y*)i — xf? + yj 
is incompressible. 
Show that any vector field of the form 

F(x, y) = fO)i + g0)j 

is irrotational. 


SURFACE AREA AND SURFACE INTEGRALS 


In Chapter 6 we were able to find the area of a surface of revolution by a single 
integral. To find the area of a smooth surface in general (Figure 13.5.1), we need a 
double integral. 


We call a function f(x, y), or a surface z = f(x, y), smooth if both partial 


derivatives of f are continuous. 


Figure 13.5.1 
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DEFINITION 


The area of a smooth surface 


z= f(x, y), (x, y)in D 


' 5 fff) Rp se 


JUSTIFICATION Let S(D,) be the area of the part of the surface with (x, y) in D,. 
S(D,) has the Addition Property, and S(D,) = 0. Consider the piece of the 
surface AS above an element of area AD (Figure 13.5.2). AS is infinitely close 
to the piece of the tangent plane above AD, which is a parallelogram with 
sides 


0 bay 
U=Axit+ —Axk, V=Ayj t+ ZAyk 
Ox Cy 


Figure 13.5.2 


The quickest way to find the area of this parallelogram is to use the vector 
product formula (Section 10.4, Problem 39), 


Area = |U x V|. 
i j k 
Oz 
Then Area = ||\Ax OQ Ax 
Ox 
Q Ay say 
yee oz. , 
= — Axi — Ax— Ayj + Ax Ayk 
Vax oy 


Teas: 
joe) 


+ (2) + 1 Ax Ay. 
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0z\? az\2 
Therefore AS = Ee + By + 1 Ax Ay (compared to Ax Ay), 


and by the Infinite Sum Theorem, 
éz\? dz\? 
= = —_ 1 dx dy. 
s= {J JE + Gl oes 
D 


EXAMPLE 1 Find the area of the triangle cut from the plane 2x + 3y+z=1 
by the coordinate planes. 


Step 1 Sketch the region as in Figure 13.5.3. 


D 
y x 

Figure 13.5.3 x 
Step 2 The plane intersects the (x, y) plane on the line 

2x + 3y = 1, yo 

Thus D is the region 

O0<x<f, o<yst5> 

Step 3 On the surface, 
z=1-—-2x — 3y, 2 U3 


az\? faz\2_— 
Then S$ = ive + (= + ldxdy 
Ox oy 
D 
= [v4 +9+1dxdy = Jia ff day 
D D 
1/2 p(l~2x)/3 /2,_9 14 
= Jia[ [ dy dx = /14 [ ee 
0 0 re) 


3 


EXAMPLE 2 Find the area of the portion of the hyperbolic paraboloid z = x* — y? 
which is inside the cylinder x? + y*? = 1. 
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Step 7 Sketch the region (Figure 13.5.4). 


Figure 13.5.4 


Step 2  D is the region 
-l<x<1l, -V/1-x*?<y</1- x? 


or in polar coordinates, 


Oz = 
oe 


2 2 
Then s=|J ued + us + 1dxdy 
Ox oy 
D 


ai /4x? + dy? + 1 dx dy. 
D 


Step 3 2x, —=-—2y 


It is easier to use polar coordinates, where 


J4x? + 4y? +1 = /4r? 41. 
2n 1 

s= | | Jar + 1dr dd. 
10) 0 


Put u=4?41, du = 8rdr, 
2n 54 
s=| [pv du ao 
o 418 


2n 1 I 
= 505°? — 140 = 2082 =) 
0) 


The line integral has an analogue for surfaces called the surface integral. 
The form of the line integral which is most easily generalized to surfaces is the vector 
form 


| F «Nas = { —Qdx+Pdy 
c c 


where N is the unit normal vector of C. This is convenient because surfaces also have 
unit normal vectors. 
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Before stating the definition we motivate it with a fluid flow interpretation, 
Remember that in the plane the line integral 


| F-Nds 


¥C 


is equal to the flux, or net rate of fluid flow across the curve C in the direction of the 
normal vector N, 
Consider a fluid flow field 


F(x, y,z) = Pi+ Qj + Rk 


and a surface S in space. Call one side of S positive and the other side negative, and 
at each point of S let N be the unit normal vector on the positive side of S. The surface 


integral 
ff F-NdS 
Ss 


will be the flux, or net rate of fluid flow across the surface S from the negative to the 
positive side (Figure 13.5.5). 

With this interpretation in mind we shall define the surface integral and then 
justify the definition. First we need the notion of an oriented surface. 


DEFINITION 


An oriented surface S is a smooth surface 


z = g(x,y) 


over a plane region D with a piecewise smooth boundary, together with an 
orientation that designates one side of the surface as positive and the other 
side as negative. (See Figure 13.5.6.) 


Zz 
vA 7 = 
¥ 
y 
x 
x An oriented surface 


Figure 13.5.5 Figure 13.5.6 


DEFINITION 


Let S be an oriented surface z = g(x, y) over D and let 


F(x, y,z) = Pi + Oj + Rk 
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be a vector field defined on S. The surface integral of F over S is defined by 


F.-NdS= + pe of + RdA, 
Hl Lee 


+ if the top ue of S is positive, — if the top side of S is negative. 


Thus a change in orientation of S changes the sign of the surface integral. 


JUSTIFICATION We show that this definition corresponds to the intuitive concept 
of flux, or net rate of fluid flow, across a surface. Suppose S is oriented so 
that the top surface of S$ is positive. 


Let B(D) be the flux across the part of S over a region D. Consider an element 
of area AD and let AS be the area of S over AD. Then AS is almost a piece of 
the tangent plane. The component of fluid flow perpendicular to AS is 
given by the scalar product F + N where N is the unit normal vector on the 
top side of AS (Figure 13.5.7). Thus the flux across AS is 


AB x F-NAS (compared to AA). 


Figure 13.5.7 


This suggests the surface integral notation 


[[F-nas 
5s 


Let us find F, N, and AS. The vector F at (x, y, z) is 
(1) F(x, y, 2) = Pi + Qj + Rk. 


From Section 13.1, one norma! vector at (x, y, z) is 


The unit normal vector N on the top side of AS has positive k component 
and length one, so 


830 
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(2) N= 


From our study of surface areas, 


AND ‘,) 2 
(3) as= / (2) + (2 +1AA (compared to AA). 
Ox ) 


When we substitute Equations 1-3 into F+ N AS, the radicals cancel out and 
we have 
Oz 


a 


= pel + R| AA (compared to AA). 
Ox oy 


Using the Infinite Sum Theorem we get the surface integral formula 


Oz oz 
B(D) = ae ee, + RdA, 
D 


EXAMPLE 3 Evaluate the surface integral 


{ F+NdS, 
Ss 


where S is the surface z = e** over the region D given by 
OSxS1 xs yp<1, 

S is oriented with the top side positive, and 
F(x, y,z) = 21+ j + 27k. 


The region is sketched in Figure 13.5.8. The first step is to find 0z/@x and 
Oz/dy. 


Figure 13.5.8 


By definition of surface integral, 
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[fronas= [f 28 -$ + 2a 


D 


1 pl 
-| | — e+ e?* "29 dy dx 
0 x 


- | —44 e771 — te 2 dx =4—e° 1 4+he-?. 
8) 


The same surface integral with S$ oriented with the top side negative has 
minus the above value. 


PROBLEMS FOR SECTION 13.5 


1 


a & ww WY 


16 


Find the area of the triangle cut from the plane x + 2y + 4z = 10 by the coordinate 
planes. 

Find the area cut from the plane 2x + 4y + z = 0 by the cylinder x? + y? = 1. 

Find the area of the surface of the paraboloid z = x? + y? below the plane z = 1. 
Find the area of the surface of the cone z = ./x? + y” below the plane z = 2. 

Find the surface area of the part of the sphere x? + y? + z* = a? which lies in the 
first octant; i.e, x = 0, y = 0,z = 0. 

Find the surface area of the part of the sphere x? + y” + z? = a* which is above the 
circle x7 + y* < b? (b < a). 

Find the surface area cut from the hyperboloid z = x? — y? by thecylinder x? + y? = a?. 
Find the area cut from the surface z = xy by the cylinder x? + y? = a?. 

Find the surface area of the part of the sphere r? + z? = a? above the circle r = acos 9. 
Find the surface area of the part of the cone z = cr above the circle r = acos 0. 

Find the area of the part of the plane z = ax + by + c over a region D of area A. 
Find the surface area of the part of the cone z = c./x? + y? over a region D of area A. 
Ene the ne area of the part of the cylinder x? wt z? = a? cut out by the cylinder 
x? + y? <a’. 

Find the surface area of the part of the cylinder x? + z? = a* above and below the 
square -b<x<b,-b<y<b(b<a). 


Evaluate the surface integral 
| fei — 3j + 4k)-NdsS, 
Ss 
where S is the surface z = x? + y?, -1 <x <1, -1 < y <1, oriented with the top 


side positive. 
Evaluate the surface integral 


[ | ei + vi + 31) -Nas 
§ 


where S is the surface z = 3x — Sy over the rectangle 1 < x < 2,0 < y S 2, oriented 
with the top side positive. 
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17 Evaluate the surface integral 


f oxi + yi — 2k)-NdsS 
s 


2 


where S is the surface z = 1 — x? — 17, x? + }* < 1, oriented with the top side positive. 


18 Evaluate the surface integral 


[foi + y2zj + zxk)>NdS 

Ss 
where S is the surface z =x + y?+2,0<5x<1,x <3} <1, oriented with the top 
side positive. 


19 Evaluate the surface integral 
[ {cect e} + zk)«NdS 
Ss 


where S is the surface z= xy, OS x <1, —x Sy Sx, oriented with the top side 
positive. 


20 Evaluate the surface integral 


| [x2 + yzj + zk 
s° 
where S is the surface z = Ja? = x? — y?, x? + y? <b, oriented with the top side 
positive (b < a). 
O 21 Show that if S is a horizontal surface z = c over a region D, oriented with the top side 
positive, then the surface integral over S is 


i a 


| [ees 2 + be, 9,25 + Rls 9, 29K) “NS = f] Rix, y, 0) dA. 


Ss D 
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Both Stokes’ Theorem and Gauss’ Theorem are three-dimensional generalizations 
of Green’s Theorem. To state these theorems we need the notions of curl and 
divergence in three dimensions. The cur] of a vector field in the plane is a scalar field, 
while the curl ofa vector field in space is another vector field. However, the divergence 
in both cases ts scalar. 


DEFINITION 


Given a vector field F(x, y,z) = Pi + Qj + Rk 


in space. The curl of F is the new vector field 


Cz ex x ey 


A A ap AR 
curl F = ( ee (¢ eli 
cy =z 


20 oP 
aan \x 
C 


This can be remembered by writing the curl as a “determinant” 
i k 

7) é 

curlF =/ =~ — = 

6x oy dz 

R 


|P Q 
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The divergence of F is the real valued function 


OP Oo é 
divF = — + g . 


ox dy 6z° 


EXAMPLE 1 Find the curl and divergence of the vector field 


F(x, y, z) = xyi + y2zj + zxk. 


ij k 
Cnn?) 
curl F = Ax oy az 
xy yZ 2x 
_ (222) _ a2), , (Ay) _ Ae), , (Alv2) _ dom)), 
oy oz Oz Ox Ox oy 
= —yi — zj — xk. 
ig OO OO OOo ge x: 
ox oy 0z 


Two interesting identities are given in the next theorem. 
THEOREM 1 


Assume the function f(x, y, z) and vector field F(x, y, z) have continuous second 
partials. Then 


curl(grad f) = 0, div(curl F) = 0. 


PROOF We use the equality of mixed partials. 


_o. of. e 
grad f Axi + ay + ao 
i j k 
6 @ @ 
curl (grad f) = ax ay a: 
ieee lie E 
Ox Ody dz 
62 Q? 2 e2 2 a2 
(Fe Ney (OE OF Ot), 
Oy dz dzdy 0z0x 0x 0z Ox dy oy dx 


The other proof is similar and is left as a problem. 


Stokes’ Theorem relates a surface integral over S to a line integral over the 
boundary of S. It corresponds to Green’s Theorem in the form 


F-Tds = [{ curl F da. 
oD fi 
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Let S be an oriented surface over a region D. The boundary of S, dS, is the 
simple closed space curve whose direction depends on the orientation of S as shown 


in Figure 13.6.1. 
The notation 


fF -Tds or p Pax + Ody + Raz 
cs 


denotes the line integral around 0S in the direction determined by the orientation of S. 


as 


x x 
Figure 13.6.1 The Boundary of S 


STOKES’ THEOREM 


Given a vector field F(x, y, z) on an oriented surface S, 


fa.F-Tas = | [ oul F Nas. 
Ss 


(See Figure 13.6.2.) 


Figure 13.6.2 


To put this equation in scalar form, let 


F = Pi + Qj + Rk, curlF = Hi + Lj + Mk. 


Then F-Tds = Pdx + Qdy + Raz, 
and if S is oriented with the top side positive, 
curlF-NdS ={-Ho —~L@ 4+ mM] da. 


Ox oy 
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Thus Stokes’ Theorem has the scalar form 


p Pdx + Qdy + Raz =[ [ oo ee 
as - é y 


x 0 


Stokes’ Theorem has two corollaries which are analogous to the Path 
Independence Theorem. 


‘COROLLARY 1 


If f(x, y,2z) has continuous second partials, then the line integral of grad f 
around the boundary of any oriented surface is zero, 


(grad fT ds = 0. 
(See Figure 13.6.3.) 


PROOF curl(grad f) = 0, so 
f grad fT ds = | { curigrad /)-N as = { fods =o. 
6s 
AY AY 


COROLLARY 2 


The surface integral of curlF over an oriented surface depends only on the 
boundary of the surface. That is, if 0S, = 0S, then 


{ feuF -N, dS, =[f curl F -N, dS. 
Si S2 


(See Figure 13.6.4.) 


grad f 


4 $s grad {+T ds =0 
Figure 13.6.3 Figure 13.6.4 


PROOF By Stokes’ Theorem, both surface integrals are equal to the line integral 


F-Tas = ¢ F-T ds. 


aS, Jas, 
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For fluid flows, Stokes’ Theorem states that the circulation of fluid around 
the boundary of an oriented surface S is equal to the surface integral of the curl 
over S. 

We shall not prove Stokes’ Theorem, but will illustrate it in the following 
examples. 


EXAMPLE 2 Let S, be the portion of the plane 
z=2x+2y—-1 
and §, the portion of the paraboloid 
z=x? 4 y? 


bounded by the curve where the plane and paraboloid intersect. Orient both 
surfaces with the top side positive, so they have the same boundary 


C = 6S, = 0S). 
Let F(x, y,z) = zi + xj + yk. 


Evaluate the integrals 


(a) jp omeen, dS,. 
St 


(b) [ few .N, ds. 
Ss 


2 


(c) br -T ds. 


By Stokes’ Theorem, all three answers are equal, but we compute them 
separately as a check. 


The regions are drawn in Figure 13.6.5. First we find the plane region D 
over which S, and S, are defined. The two surfaces intersect at 


2x + 2y—1 =x? +4 y?’, 
(x -—1P +(y-1P =1. 


Figure 13.6.5 


13.6 THEOREMS OF STOKES AND GAUSS 


So D is the unit circle with center at (1, 1) shown in Figure 13.6.6; that is, 


O<x <2, 1-./1-—(x-1P <y<1+/1-—(@-— 1). 


y 
x 
Figure 13.6.6 
Next we compute curl F. 
j &k 
OO. 8 
ewlF =|— — ~—]=i+jt+k. 
ex Gy éz 
Zo oR 
(a) On the surface z = 2x + 2y — 1, 
Oz oz 
—=2, — =2. 
Ox oy 


Thus | cone, as, = {[-2-2+ 1dA= -3f [aa = —3n 
D D 


Si 


(b) On the surface z = x? + y’, 


Thus (eee 1dA 
So 


1+/1—(x-1)? 
-[f — 2x — 2y + ldydx 
1-/1-(x-1)? 


aieneae 
-{ —2xy — y? +y| dx 
0 1-J/1-@— 12 
2 secs, 
al dy. /1 — Ge — 1)? ~ 2/1 — (& = 1)? dx 
0 


= —3n. 
(c) The boundary curve C = 0S, = 0S, is a space curve on the plane 
= = 2x + 2y — 1 and over the circle 

(x — 1)? +(y—- 17 =1. 

Thus C has the parametric equations 
x=1+cos0, y=1+sin0, z=2cos0+2sin0+ 3, 0<0<2z. 

Then dx = —sin @ dé, dy = cos dé, 

dz = (—2sin@ + 2cos 0) dé. 
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pF Tds =} zdx + xdy + yde 
c © 


- [tacos + 2sin 0 + 3)(—sin 0) + (1 + cos 0) cos 0 
° + (1 + sin 0)(—2sin 0 + 2 cos 0)] dO 
= [va + 3cos@— S5sin0 — Ssin? 0)d@ = —3z. 
Notice that (a) ee much easier than (b) or (c). 


Gauss’ Theorem shows a relationship between a triple integral over a region 
E in space and a surface integral over the boundary of E. It corresponds to Green’s 


Theorem in the form 
) F-Nds = divF dA. 
eD 


D 


Before stating Gauss’ Theorem, we must explain what is meant by the surface 
integral over the boundary of a solid region E. In general, the boundary of E is made 
up of six surfaces corresponding to the six faces of a cube (Figure 13.6.7). Sometimes 
one or more faces will degenerate to a line or a point. 


Figure 13.6.7 


The top and bottom faces of E are (x, y) surfaces, that is, they are given by equations 
z = c(x, y). However, the left and right faces of E are (x, z) surfaces y = b(x, z), while 
the front and back faces of E are (y, z) surfaces of the form x = a(y, z). Surface integrals 
over oriented (x, z) and (y, z) surfaces are defined exactly as for (x, y) surfaces except 
that the variables are interchanged. 

In the following discussion E is a solid region all of whose faces are smooth 
surfaces. 


DEFINITION 


The boundary of E, GE, is the union of the six faces of E oriented so that 
the outside surfaces are positive. The surface integral of a vector field F(x, y, z) 


over OE, 
[{F-Nas, 
CE 


is the sum of the surface integrals of F over the six faces of E. (See Figure 
13.6.8.) 
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Boundary of E 


Figure 13.6.8 


We are now ready to state Gauss’ Theorem. 


GAUSS’ THEOREM 


Given a vector field F(x, y, z) and a solid region E, 


[[p-nas= | [aivrar, 
CE E 
This equation may also be written in the form 


aP 00. OR 
[[renas=[[ [S++ 5 av 


6E E 


Gauss’ Theorem is sometimes called the Divergence Theorem. 

For fluid flow, Gauss’ Theorem states that the outward rate of flow across 
the boundary of E is equal to the integral of the divergence over E (Figure 13.6.9), 
As in the two-dimensional case, the divergence is the rate at which the density is 
decreasing. 


Figure 13.6.9 


The following corollary is another analogue of the Path Independence 
Theorem. 
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COROLLARY 3 


If F(x, y,2z) has continuous second partials, the surface integral of curlF 
over the boundary of E is zero. In symbols, 


[ [ cut -Nas =) 
CE 


PROOF Since div(curl F) = 0, 
lf curlF -NdS = [ ff aiveunn av = | ff OdV = 0. 
E he E 


EXAMPLE 3 Use Gauss’ Theorem to evaluate the surface integral 


( } F-NdsS, 

oe 
where F(x, y,z) = e*i + ej + xyzk 
and E is the unit cube in Figure 13.6.10. 


O<xsil, O<sy<l, O<sz<l. 


Figure 13.6.10 


By Gauss’ Theorem, 


[yeeses || aed 
=|{fete+xyav 


E 
1 i 1 
=| | [ e+e + xydzdydx 
0 v0 ¥0 
1 1 
at i e+e + xy dy dx 
6 40 
i 
=| e+e—1+4xdx 


0 
= 1 
= 2e — 5. 
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PROBLEMS FOR SECTION 13.6 

In Problems 1-6, find the curl and divergence of the vector field. 
1 F(x, y,z) = x71 + y?j + 27k 

F(x, y, z) = x cos zi + ysin zj + zk 

F(x, y, 2) = (x + y + Zi + (y + 2)jj + zk 

F(x, y, z) = yzi + xzj + xyk 

F(x, y, z) = xe'*7i + ye**4j + ze*k 

F(x, y,z) = yit+ xj +k 


AA Un & WS LV 


div(curl F) = 0. 
8 Given a function f(x, y, z) with continuous second partials, show that 
af af af 
ax? by oy? + dz?" 


div(grad f) = 


9 Use Stokes’ Theorem to evaluate the surface integral {/; curlF -NdS where S is the 
portion of the paraboloid z = 1 — x? — y? above the (x,y) plane and F(x, y, z) = 


xy*i — x”yj + xyzk. (S is oriented with the top side positive.) 
10 Use Stokes’ Theorem to evaluate the line integral 


) (Qi + zj — xk) - T ds 
és 


where S is the portion of the plane z = 2x + 5y inside the cylinder x? + y? = 1 oriented 


with the top side positive. 
il Use Stokes’ Theorem to evaluate the line integral 


) (ax + by + cz)\it+j+k)-Tds 
as 


where S is the portion of the plane z = px + gy + r over a region D of area A, oriented 


with the top side positive. 
12 Use Stokes’ Theorem to show that the line integral 


> (P@)i + Q0)) + R(z)k) + T ds = 0 
6S 


for any oriented surface S. 


13 Use Gauss’ Theorem to compute the surface integral 
[fori + y*j + 27k)-NdS 
@E 


where E is the rectangular boxO <x<a4,0Sy<b,0<zc. 


14 Use Gauss’ Theorem to compute the surface integral 


{f (2xyi + 3xyj + ze**k)-NdS 
GE 


where E is the rectangular boxO <x<1,0Sysl10<z<l. 


15 Use Gauss’ Theorem to evaluate 
{ (xi + 2yj + 3zk)-NdS 
€E 


where E is the regionO<x<10Sysx,0SzsSxty. 


Prove that for every vector field ‘F(x, y, z) with continuous second partials, 
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16 Use Gauss’ Theorem to evaluate 
[ fos + y°9j + 22k)-NdS 
CE 
where E is the sphere x? + y? + 2? <4. 


17 Use Gauss’ Theorem to evaluate 
[[ wv +y+2)(i+j+k-Nds 
cE 
where E is the hemisphere 0 < z < ./1 — x? — y?. 


18 Use Gauss’ Theorem to evaluate 
J (xy?i -+ yzj + x?zk)«NdS 
cE 
where S is the cylinder x? + yp? <1,0<2<4. 


19 Use Gauss’ Theorem to evaluate 
[ ce cos? zit ysin? zj + ./x? + y? zk) -NdS 
CE 
where E is the part of the cone z = | — \/x? + y* above the (x, y) plane. 


EXTRA PROBLEMS FOR CHAPTER 13 


I Find the derivative of z = cos x + sin y in the direction of the unit vector U = cos ai 
+ sin aj. 
2 Find grad f and fy if 
{(% y) =coshxsinhy, Us= = zs J 
we) 


3 Find grad f and fy if 
S(x,y) = e, U = cos ai + sin aj. 
4 Find the derivative of z = In (x? + y?) at the point (—1, 1) in the direction of the unit 
vector U = ai + bj. 
Find a unit vector normal to the surface z = xy at the point (2, 3, 6). 


Evaluate the line integral 
{ (cos xi — sin yj) dS 
c 


where C is the curvex =?,y=P,0<t< 1. 


Dy god 
[ilo tar] ® 


where C is the rectangular curve from (1, 2) to (4, 2) to (4, 4). 


7 Evaluate the line integral 


8 Evaluate the line integral 
{ (xi + yj + zk)-dS 
c 
where C is the line x = 2t,y = 3f,z2 = -4,0S¢°S 1. 
9 Find the work done by the force F = yi + x?j acting once counterclockwise around 


the circle x? + y? = 1. 


10 Find a potential function for y cosh xi + sinh xj. 


11 
12 


13 
14 


15 


16 


17 


18 
19 


20 


21 


22 
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EXTRA PROBLEMS FOR CHAPTER 13 


Find a potential function for (y In y + In x)i + (xIny + x)j. 
Solve the differential equation 

(2x — 6x?y + y)dx + (—2x3 + 3xy? + 1)dy =0. 
Solve the differential equation e~” sin x dx + (e~* cos x + 3y) dy. 


Use Green’s Theorem to evaluate the line integral 
$ sin x sin ydx + cos x cos y dy, 
aD 


D:t/6<x<2/3,7/6 Sy < 2/3. 
Use Green’s Theorem to evaluate the line integral 
2xy? dx + 3x*y3 dy, D:0Sx<1,x? < y< 2x. 
6D 
Use Green’s Theorem to find the area of the region bounded by the parametric curve 
x = acos? 6, y = asin® 0, 0<6<2z. 


Find the area of the part of the surface z = x” + y which lies over the triangular region 
O<xsLlOsy<x. 


Find the area of the part of the surface z = xy which is inside the cylinder x? + y? = 4. 


Evaluate the surface integral 
[fc + yj + 2k)-NdS, 
AY 


where S is the upper half of the sphere x? + y? + z? = 1, oriented with the top side 
positive. 
Find the curl and divergence of the vector field 
F(x, y) = xe’i + ye*j. 
Find the curl and divergence of the vector field 
F(x, y, z) = xyzi + xy?z3j + x?yzk. 
Use Gauss’ Theorem to evaluate the surface integral 
{{ (xy?i + yz?j + x?yk)-NdS 
OE 
where E is the region x? + y?><1,x7+y?*<z<1. 
The gravitational force of a point mass m, acting on another point mass m, has the 
direction of the vector D from m, to m, and has magnitude proportional to the inverse 
square of the distance |D|. Thus 
_ emym,D 
(DP 


where c is constant. Use the Infinite Sum Theorem to show that the gravitational force 
of an object with density h(x, y, z) in a region E on a point mass m at (a, b, c) is 


F = Pi + Qj + Rk, 


cmh(x, y, z)(x — a) 

where P= MN [« — a? +(y— bP + 2 — oP? dV, 
emh(x, y> z)(y <= b) 

g - JJ (x — a2 + (y — bY +e - path 


cmh(x, y, z)(z — ©) 
R= Ne ae +(y_bP + opr” 
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Suppose z = f(x, y) is differentiable at (a,b). Prove that the directional derivatives 
Sy(a, b) exist for all U. (See also extra Problem 36 in Chapter 11.) 


Let U =cosai + sinaj. Suppose that z = f(x,y) has contmuous second partial 
derivatives. Prove that the second directional derivative of f in the direction U is given 
by 


; i ae esr : af, 
x, y) = 5 c0s*a +2 cos «sing + =“5sin* a. 
foul ¥) x2 + Ox dy ay? 
Second Derivative Test for two variables. Suppose 
(a) (x,y) has an interior critical point (a, b) in a rectangle D. 
agp ape can 
: Cf ef écf . 
b) Throughout D, -;,—5,=— are continuous and 
(6) 8 éx?* @y?? Ax ay 
az adr arp a2 q2y\2 
¢ é Cee fe 
eS 5 9, ea mist ~ (al > 0 
ex y ex? ey Oxdy 


Prove that f has a minimum in D at (a,b). Hint: Use the preceding problem to show 
that all the second directional derivatives fyy(x, ») are positive so that the surface 
z= f(x,y) has a minimum in every direction at (a,b). In the case cosxsina > 0, 
use the inequality 


af a ee 

Os [5 z7cos%~— / x5 sin 4 
N Cx Nicy 

and use a similar inequality when cos x sin x < 0. 


Given a sphere of mass m, and constant density, and a point mass m, outside the 
sphere at distance D from the center. Show that the gravitational force on m, is the same 
as it would be if all the mass of the sphere were concentrated at the center. That is, 
F points toward the center and has magnitude 
cmym, 

D? 
Hint: For simplicity let the center of the sphere be at the origin and let m, be at the 
point (0, 0, D) on the z-axis. Let the sphere have radius b and density h, so 


IF| = 


h = m,/volume = 3m,/4nb°, ob < D. 


By symmetry the i and j components of the force are zero. Use spherical coordinates 
to find the k component, 


_ em ah+(z — D) 
: i Weescescet auald 


2x pb px cmh(p cos @ — D)p? sing 
7 i. ie I, [p? + D? — 2Dpcos 6}? 
A region D in the plane has a piecewise smooth boundary 0D and area A. Use Green’s 
Theorem to show that an object with constant density k in D has center of mass 
wees x? dy, p= aah y? dx. 
2A Jap : 2A Jap 
Show that the object in the preceding exercise has moment of inertia about the origin 


k 
T= — yrdx + x3 dy. 
3 Jop 


do dp de. 


Use the Infinite Sum Theorem to show that the mass of a film of density p(x, ») per unit 
area on a surface z = f(x, y), (x, py) in D, is 


n[f 


+1 p(x, y)dx dy. 
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EXTRA PROBLEMS FOR CHAPTER 13 


Show that the volume of a region E is equal to the surface integral 


V= $f fos + yj + zk)-NdsS. 
GE 
Show that the gravity force field of a mass m at the origin, 
m xi+ yj + zk 
F = , 
ONDE EE aby yea 


is irrotational (except at the origin). Use Stokes’ Theorem to show that 


F(, y,z)+Tds =0 


és 
where S is any oriented surface not containing the origin. 


Show that for any smooth closed curve C around the origin, 


-y * 
dx + dy = 2n. 
bare epg te 


Assume for simplicity that C has the parametric equation 
C:r=f(0#),0<0< 2n where 0 < f(0), f(0) = f(2z). 
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